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Abstract 



Effects of Coulomb interaction on persistent currents in disordered one- 
dimensional rings are numerically investigated. First of all effectiveness of the 
Hartree-Fock approximation is established on small systems. Then the calcu- 
lations are done for systems with 40 electrons in 100 sites. It is found that 
the amplitude of the average persistent current in the diffusive regime is sup- 
pressed as the strength of the Coulomb interaction increases. The suppression 
of the current is stronger in larger rings than in smaller ones. The enhance- 
ment of the current by the electron-electron interaction was not observed in 
the diffusive regime. 
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About ten years ago, Biittiker et al. [|T| suggested that equilibrium persistent current 
could exist in a mesoscopic normal metal ring pierced by a magnetic flux, if the size of 
the ring is so small that the coherence of electrons are kept in the whole system. Later the 
existence of the persistent current was conflrmed by three experiments. The flrst experiment 
was done on many isolated copper rings, the second one |^ on a single gold ring, the 
third one on a GaAs/AlGaAs semiconductor. In the third experiment the system is in 
the ballistic regime, and the size of the persistent current is in reasonable agreement with 
theories: The current is of the order of Jq = ev-p/L, where v-p is the Fermi velocity and 
L is the circumference of the ring. On the other hand, in the former two experiments the 
system is in the diffusive regime. In this case simple theories which neglect the electron- 
electron interaction had predicted the persistent current of the order of / ~ Iq/M where 
M is the number of transverse channels However, the persistent current observed 

was greater than that. Especially in the second experiment the magnitude of the current 
was comparable to the clean case: / ~ Jq. We might think that the discrepancy comes 
from inadequacy of the approximations employed. However, numerical calculations done by 
Montambaux et al. |^,^ for non-interacting electrons in disordered rings also cannot explain 
the large magnitude of the current as seen in the second experiment. Thus it is natural to 
consider that the discrepancy between theoretical calculations and experimental data is not 
due to the approximation employed in the theory but due to a defect of the model where 
electron-electron interaction is neglected. 

Models with electron-electron interaction have been investigated both analytically and 



numerically. In some analytic calculations |10-12| the diagrams called cooperon which most 



contribute to the flux dependence of the energy were considered and resulted in larger 
current than that in non-interacting case. However, it was found that the effect of higher 
diagrams suppress the current [Q. In numerical calculations effects of the electron-electron 
interaction were investigated by exact diagonalization of its Hamiltonian for small size rings. 
In the case of a one-dimensional discrete-lattice ring the interaction, which are both long 



ranged [0 and short ranged ||15[, is found to suppress the average persistent current. On the 



contrary, starting from continuum model it was suggested that the average persistent 



current was enhanced by the electron-electron interaction. Thus, in both analytical and 
numerical calculations the role of the electron-electron interaction for the persistent current 
is controversial. 

The difficulty in the present problem is that there is no guiding principle for the analytical 
calculation. It is not clear which diagrams we should take into account. On the other hand 
numerical investigations, where we can obtain results free from approximations, are limited 
to very small systems. Thus it is desirable to do numerical calculation for much larger 
systems. This is what we seek in this paper: to numerically investigate the effects of the 
Coulomb interaction on the persistent current in a larger sample than ever investigated and 
to find whether the results will be different from those of the small size or not. In order to 
study larger systems we use the Hartree-Fock approximation (HFA). We first compare the 
results of the HFA with those by the exact diagonalization for small systems. We find that 
the HFA gives qualitatively correct results. So we apply HFA to systems with 100 sites for 
which exact diagonalization is impossible. 

We consider a one-dimensional disordered lattice ring with A^^ sites pierced by a magnetic 
fiux in which electrons mutually interact with long-range (1/r) Coulomb repulsion. For 
simplicity we neglect the spin degrees of freedom of electrons. We adopt the tight-binding 
model for kinetic term. The impurity potential is introduced by random site-energy Si which 
has random uniform distribution with width W {—W/2 < Si < W/2). Thus, 

n = -t ^(e'^aj+itti + e"'^ajai+i) + ^ e^alai 



2^Ar,|2sin[^(z-j)]| ' ' ' ^ ^ 

where Oj (a|) is an annihilation (creation) operator of a spinless fermion at z-site. The 
Ng + 1-th site is identical to the first site. The effect of a flux (0) piercing a ring is introduced 
through the phase factor e'^ which is gained when a electron hops to the neighboring site, 
thus 9 = 2'K(f)/Ns4>o, where 0o = h/e. Coulomb repulsion is parameterized by V{> 0). 
In the HFA H is replaced by T^hf, 



i=l i=l 



1 „ V {{a^jaj)a\ai — {a},jai)a\aj) 
^2tr^ |2sin[^(^-j)]l ' 
where (■ ■ ■) denotes the expectation value with respect to self-consistent HF eigenstate l^') = 
^1^2 ■ ■ ■ ^]ve|0)' where 6„ is the electron-annihilation (creation) operator of the n-th 
eigenstate of TYhf and is the total number of electrons. The self-consistency is achieved 
by iteration. Once the Hamiltonian is solved the persistent current is calculated by the 
formula, 

- (3) 

where Eg{(j)) is the ground state state energy at flux 0. 

In the following we calculate the current at several values of W/t and V/t. We mostly 
calculate the current in the diffusive regime. In one dimension the localization length ^ is 
given by 

? = ^ {W«2n), (4) 

where a is the lattice constant [p!7H l9|. Equations above were derived under the condition 
of the half-filled-band case for the tight-binding model. The system investigated here is not 
half-filled. However, we use these equations to estimate the localization length, since it is 
not far from half-filled. We are mainly interested in the diffusive regime and consider the 
value of W such that ^ > L = Nga. 

We first examine validity of the HFA. For that purpose we calculate the persistent current 
for systems of 4 electrons in 10 sites by the two methods: the exact diagonalization and 



the HFA. Average over site energy randomness is performed over ten samples |20|. Flux 
dependence of the persistent current at W = t and for = to 20t is shown in Fig.l for 
half period, < < 0o/2; the current is an odd function of 0, and it is periodic with 
period 0o- The current is normalized by the maximum persistent current for clean {W = 0) 



and non-interacting (V = 0) system, Jq = {2et / Ngh) sm^Neii / Ng) [|T9|. In these systems 
the amphtude of the persistent current is suppressed as the Coulomb interaction becomes 



larger. The result by the exact diagonalization is in agreement with that in Ref. where 
a system with 5 electrons in 10 sites is investigated. Although the suppression is stronger, 
we see the HFA gives qualitatively similar behavior. The comparison is done at other values 
of W also. Figure 2 summarizes the results. The vertical axis shows the magnitude of the 
current / at = 0o/4 divided by /dean which is the value of the current at = 0o/4 in 
the clean, non-interacting system: /dean = /o sin(7r/2A^s)/ sin(7r/A^s) — /o/2 [l^. The plot of 
non-interacting case (circle) is almost on the curve represented by exp{—L/^{W)). The HFA 
results reproduce those of the exact diagonalization qualitatively: The persistent current is 
suppressed by the Coulomb interaction in the diffusive regime, W < 3. It is remarkable that 
the slight enhancement of the current in the localized regime {W > 3) is also reproduced. 

Now that we have found that the HFA gives qualitatively correct behavior of the per- 
sistent current, we investigate larger systems where the exact diagonalization is impossible. 
Specifically we consider systems with 40 electrons in 100 sites. Figure 3 shows the flux 
dependence of the persistent current at W = 0.25t, and Fig.4 shows W dependence of the 
current. In these figures average over 10 samples are shown. For this 100-site systems the 
sample dependence is small due to self-averaging. Therefore 10 samples are enough. These 
results also show the suppression of the current with increase inV. In this system ,^ ~ L at 
14^ ~ t, so only the diffusive regime is considered. The non-interacting case is also almost 
on the curve exp{—L/W{^)). Compared to the 10-site system it is seen that the effect of V 
is larger in the present larger system. 

From the results stated above we can conclude that the Coulomb interaction between 
electrons cause the reduction of the persistent current in the experimentally relevant diffusive 



regime. The results found in Refs. [0 and [|T5| are confirmed for larger systems. In Ref. 



it was suggested that this suppression was attributed to the Mott-Hubbard transition. 
If so, it should occur in a clean {W = 0) system, and in fact it was observed in Refs. 
1^ and ||15|. It is also seen in our results. Fig. 2 (a). On the other hand, in the HFA 



the suppression of the current does not occur in a clean sample at any filling and V, since 
the eigenstates in the case of V=0 are still the eigenstates and the ground state is also 
unchanged. Nevertheless, in a dirty system, the current is suppressed as V increases both in 
the case of the exact diagonalization and the HFA. Furthermore the Coulomb effect is larger 
in disordered systems in the diffusive regime. This seems to teach us that the Mott-Hubbard 
transition is not essential to the suppression of the current in a dirty system. 

Our results are in disagreement with those in Ref . |jl6| . That may be considered due to the 
difference in models, discrete or continuum, as some people |]2ll , ^ pointed out. However, in 
view of the fact that the system with more lattice points approaches the continuum system, 
our results that the Coulomb interaction suppress the current more in the larger system cast 
some doubt to the interpretation of the difference. Another candidate for the reason of the 



enhancement in Ref. ||T6| would be truncation of the single electron states adopted there. 

We cannot explain the experimental results by our model and calculations. It may be 
possible that the enhancement of the current will occur in multi-transverse channel systems. 
These systems cannot be investigated by the exact diagonalization because of the limitation 
in memories. However, we can treat such a system by the HFA. Results of such a calculation 
will be published in near future. 

In summary we have first established the effectiveness of the HFA for the present model. 
With this approximation it has become possible to investigate systems with more than 100 
sites even in the presence of the Coulomb interaction. We have found that the Coulomb 
interaction suppress the persistent current. 

One of the authors (D. Y.) thanks Aspen Center for Physics where part of this work was 
done. This work was supported by a Grant-in Aid for Scientific Research on Priority Area 
"Computational Physics as a New Frontier in Condensed Matter Research" (04231105) from 
the Ministry of Education, Science and Culture. 
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FIGURE CAPTIONS 



FIG.l. Current I / Iq vs flux 0/0o for a 10-site 4-electron ring dXW — t^ and V — Q (circle), 

5t (square), lOt (diamond), and 20t (triangle), where Jq = {2et / N sK) sin{N eii / N s) is the 
maximum persistent current for clean, non-interacting system. The panel (a) shows the 
results by the exact diagonahzation and panel (b) shows those by the HFA. 

FIG. 2. Current ///dean vs randomness {W) for a 10-site 4-electron ring (a) by the ex- 
act diagonahzation and (b) by the HFA for y = (circle), bt (square), lOi (diamond) and 
20i (triangle), where /dean = /o ^va.{ii /2Ns)/ sin(7r/A^s) is the value of the current at = 0o/4 
in the clean, non-interacting system. The lines are guides to the eye. 

FIG. 3. Current ///q vs flux 0/0o for a 100-site 40-electron ring in the case of the HFA 
aA,W = 0.25t, and ^ = (circle), t (square), 3t (diamond) and bt (triangle). 

FIG. 4. Current I/Idean vs randomness {W) for a 100-site 40-electron ring by the HFA 
for V = (circle), t (square), 3t (diamond) and bt (triangle). The lines are guides to the 
eye. 
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